NACA TN No. 1767 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1767 

THE APPLICATION OF GREEN'S THEOREM TO THE SOLUTION 

OF BOUNDARY -VALUE PROBLEMS IN LINEARIZED 

SUPERSONIC WING THEORY 

By Max. A. Heaslet and Harvard Lomax 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 



Washington 
April 1949 


GQI t ''" 3 — ST A’ ( !i 1 LI8!fe*rR.Y 

APR 13 1949 


BUSINESS, SCIENCE 
& TECHNOLOGY DEP’T, 






NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 1767 

THE APPLICATION OF GREEN’S THEOREM TO THE SOLUTION 
OF BOUNDARY-VALUE PROBLEMS IN LINEARIZED 
SUPERSONIC WING THEORY 1 
By Max. A. Heaslet and Harvard Lomax 


SUMMARY 

Following the introduction of the linearized partial differenr- 
tial equation for nonsteady three-dimensional compressible flow, 
general methods of solution are given for the two— and three- 
dimensional steady— state and two-dimensional unsteady— state equations . 
It is also pointed out that, in the absence of thickness effects, 
linear theory yields solutions consistent with the assumptions made 
when applied to lifting— surface problems for swept— back plan forms 
at sonic speeds. The solutions of the particular equations are 
determined in all cases by means of Green’s theorem and thus depend 
on the use of Green’s equivalent layer of sources, sinks, and 
doublets. Improper integrals in the supersonic theory are treated 
by means of Hadamard’s "finite part" technique. 

Four applications of the general solutions are given: First, 

the angle— of— attack load distribution for a supersonic, yawed, 
triangular plate with subsonic leading edges is determined. Second, 
downwash is calculated along the center line in the plane of the 
unyawed triangular wing. Third, the growth of load distribution is 
presented for subsonic and supersonic two-dimensional flat plates 
either starting from rest at a uniform velocity or experiencing an 
abrupt angle— of— attack change. The transient effects on lift— curve 
slope are then calculated. Finally, the load distribution and lift- 
curve slope of a specific swept— back lifting surface are determined 
at a free— stream Mach number of one. 


presented by Dr. Heaslet at the VII International Congress of 
Applied Mechanics, September 19^8, London. 


2 


MCA TN No. 1767 


INTRODUCTION 

If the effects of viscosity are assumed, small and shock— free 
compressible flow is considered, the velocity field about a two— or 
three-dimensional body placed in a uniform free stream is irrota— 
tional and thus possesses a velocity potential. In the determination 
of the pressures exerted on such a body or in the calculation of the 
induced velocity components, the theoretical aerodynamicist is 
concerned essentially with finding the velocity potential of the flow 
field and, thus, must determine the solution of a second— order non- 
linear partial differential equation subject to certain boundary 
conditions. The known mathematical difficulties that arise in the 
treatment of such a problem make it expedient to resort to simplify- 
ing assumptions. In applied aerodynamics, however, efficiency of 
flight at high speeds has focused attention on bodies inducing 
relatively small velocities throughout the field of flow and, as a 
consequence, the demands of engineering furnish a guide for the 
mathematical simplification of the theory. The so-called linearized 
theory of compressible flow was developed to solve such problems and, 
although considerable work of a more precise nature has been presented 
in two dimensions, a large amount of investigation in unsteady flight 
and in three-dimensional wing theory remains to be completed within 
the framework of the simplifying conditions. 

The present paper is restricted to a discussion of wing theory 
subject to the assumptions of linearized compressible flow. It 
therefore employs solutions of La.place’s equation and the wave 
equation for cases where the boundary conditions are specified in 
the plane of the wing. Attention will be directed primarily to the 
analysis of steady— state conditions although an equivalence will be 
established between the two-dimensional differential equation con- 
taining the time variable and the equation applying to three- 
dimensional supersonic wing theory. Solutions in all cases will be 
obtained through the use of Green’s theorem and the resultant concept 
of Green’s equivalent layer of sources, sinks, and doublets. The 
correspondence between the theoretical development for subsonic and 
supersonic speeds is particularly useful since experience related to 
analysis in either flight regime is more readily transferred. 

In view of the widespread use of sources, sinks, and doublets 
in low— speed studies and the fact that the earlier applications to 
supersonic wing theory by Prandtl (reference 1) and Schlichting 
(reference 2 ) corresponded to the use of Green’s equivalent layer, 
it is notable that later emphasis has shifted to other methods of 
solution. Sources alone were used by Puckett (reference 3) to create 
symmetrical nonlifting wings and were also applied to the study of 
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lifting triangular wings with supersonic leading edges, but the use 
of source, sink, and doublet sheets has not been as extensive as 
might have been expected. This anomaly is even more apparent in view 
of the vast mathematical and physical literature centering around 
the use of Green’s theorem. One possible explanation may stem from 
the fact that the interest of the mathematician and physicist in the 
wave equation has arisen in connection with problems in acoustics, 
optics, and vibrating membranes. Such problems introduce boundary 
conditions of the Cauchy type, that is, initial conditions need to 
be known both for the unknown function and its rate of change. The 
supporting surface for such boundary conditions cuts the character- 
istic cone of an arbitrary point in a closed curve and has been 
called by Hadamard (reference 4 ) a duly inclined surface. In aero- 
dynamics the supporting surface is nonduly inclined and cuts the 
characteristic surface or Mach cone along the arc of a hyperbola and, 
as a result, the problem is no longer of the Cauchy type and the 
analysis becomes similar to that used in subsonic theory in the solu- 
tion of Laplace’s equation. Prior to the interest of the theoretical 
aerodynamicist in supersonic wing theory, it appears that little 
attention in application was paid to this type of solution. 

The material presented here is divided into two main divisions: 
Analysis and Applications. In the first part of the Analysis 
division, the linearized differential equation for nonsteady com- 
pressible flow is given together with the underlying assumptions made. 
Specific forms of this equation for two— and three-dimensional steady 
states and two-dimensional unsteady states are then considered. It 
is pointed out in particular that for swept— back lifting surfaces 
linearized theory yields consistent solutions at a free— stream Mach 
number of one although the analysis of arbitrary thickness distribu- 
tions is not possible. Following the various equations, Green’s 
theorem is applied to find, in terms of the known boundary conditions, 
the desired solution by means of source and doublet distributions. 

Applications of the general methods are confined to four problems . 
As an example of the manner in which angle— of— attack load distributions 
are determined for a lifting flat plate, the case of a yawed triangular 
wing with subsonic leading edges is solved. Doublet distributions are 
then applied in the second problem to the calculation of downwash 
behind the same wing in an unyawed position. Third, the growth of 
load distribution with time is derived for a supersonic two-dimensional 
flat plate either experiencing a sudden sinking motion or starting 
from rest at a uniform velocity. Such distributions are of value in 
the calculation of indicial lift functions and can be used, together 
with Duhamel’s integral, in the study of certain dynamic maneuvers. 

The final application considers at a Mach number equal to one the case 
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of a swept— back lifting surface with tips normal to the free— stream 
direction. 

All important symbolo are defined in the appendix. 


ANALYSIS 

The Partial Differential Equations 

Basic differential equation .— Consider an aerodynamic body 
flying at an arbitrary Mach number Mq in air initially at rest. 

If a Cartesian coordinate system x, y, z is fixed relative to the 
body, the body may then be assumed stationary and situated in a free 
stream with the same Mach number. If the free— stream velocity vector 
is parallel to and in the direction of the position x axis and if cp 
denotes the perturbation velocity potential for isentropic flow, the 
linearized partial differential equation for cp may be written in 
the form 

(Mo 2 — l) c Pxx ~ Tyy — 9zz + ~ ~ 9tt + Txt = 0 (l) 

a o a o 

where a Q is the velocity of sound in the free stream and t denotes 
time. 


The assumptions underlying the derivation of equation (l) have 
been stated in numerous places but are not always obviously compatible. 

It is assumed here that the ratios are small compared to 

, v o v 0 v 0 

one, where u, v, w are induced velocity components and V Q is the 
velocity of the free stream; moreover, 


7-1 

2 


Mo 2 



u 2 +v 2 +w 2 

~ iT 2 

Y o 


«1 


and, finally, the velocity gradients at a given point of the flow 
field are all of similar magnitude. 

Special cases .— The particular forms of equation (l) to be 
considered are given in table I. In the steady-state equations the 
original independent variables are retained; the two-dimensional 
unsteady-state equation has, however, new variables defined by the 
relations 
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x* = x - aoMot, z* = z, t» = a 0 t (2) 


and consequently the boundary conditions for any particular example 
will be subject to the same transformation. In all the equations 
the constraints imposed by the linearization permit, for problems in 
wing theory, the boundary conditions to be specified in the plane of 
the wing. This plane shall arbitrarily be taken to be z=0. 


TABLE I.- LINEARIZED PARTIAL DIFFERENTIAL EQUATIONS OF 

COMPRESSIBLE FLOW 



Steady State 



Two J 

' (1 -Mq 2 ) <T>xx + CP ZZ = 0 , 

Mq<1, 

(A) 

Dimensions 1 

>- (Mo 2- 1) 'Txx — 9zz = 0, 

Mo > 1, 

(B) 


f (i-Mo 2 ) + Vyy + <Pzz 

= 0, M q <1, 

(c) 

Three -< 

1 (Ho 2- 1) Txx — c Pyy — ^zz 

= 0 , Mq > 1 , 

(D) 

Dimensions 

L'Pyy + ^zz = 0 

M 0 = 1, 

(E) 


Unsteady State 



Two 

Dimensions 

t*t» - Tx» x » ~ T'z*z* = 

0 

(F) 


The Mach number range for which the equations are valid cannot 
be prescribed a priori since induced velocities are functions of wing 
geometry and angle of attack. We can say, however, that for certain 
configurations at small angles of attack the equations and their 
solutions are consistent with the assumptions. In particular, three- 
dimensional lifting surfaces with sufficient sweepback yield solutions 
of this class at Mq= 1. The differential equation shows that in this 
case the boundary conditions need only be specified along strips in 
the transverse direction. The surfaces of the Mach cones also are 
normal to the free-stream direction so that any disturbance point 
makes itself felt at all points not upstream of it. Since for these 
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lifting surfaces the disturbances do not become excessive at Mq= 1, 
we have a specific kind of lateral strip theory that yields formal 
solutions compatible with the assumptions made. 


Boundary Conditions 

Steady state .— The boundary conditions are given in the z=0 
plane and in the case of two-dimensional theory the wing is assumed 
to extend infinitely, parallel to the y axis. As a convenience in 
notation, two subscripts will be introduced: the first, u, denotes 
conditions on the upper surface of the wing, that is, the limit of 
the function as z approaches zero through positive values; the 
second, 2, denotes conditions on the lower surface of the wing, that 
is, the limit of the function as z approaches zero through nega- 
tive values. 

Four types of boundary conditions arise in practice: 

1. Symmetrical nonlifting wing (boundary— value problem of the 

first kind) .- The conditions w u = W7 = 0 hold over all 
of the xy plane except for the region occupied by the wing. 
On the wing, the relations 2w u = — 2w 1 = Aw 0 = f(x,y) 
are given, the function f(x,y) being determined by the 
geometry of the configuration. Over all of the xy plane, 
Au 0 = u u — u^ =0 applies. 

2. Lifting surface with specified loading (boundary value 

problem of the first kind) .— The condition 
Au 0 = u u — ui = 0 holds over the xy plane except for the 
region occupied by the wing. On the wing, the relations 
2u u = - 2u 1 = Au 0 = f(x,y) are given, the function f(x,y) 
being determined by the specified loading. Over all of 
the xy plane, Aw 0 = 0 applies. 

3. Lifting surface with specified camber and angle of attack 

(boundary— value problem of the second kind) .— The 
condition Au 0 = 0 holds over the xy plane except for 
the region occupied by the wing. On the wing the 
relation w = f(x,y) is given, the function f(x,y) 
being determined by the given camber, twist, and angle 
of incidence. Over all of the xy plane* Aw 0 = 0 applies. 

4. Symmetrical wing with specified pressure distribution 

(boundary— value problem of the second kind) .— The 
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condition Aw 0 = 0 holds over all of the xy plane except 
for the region occupied by the wing. On the wing the 
relation u = f(x,y) is given, the function f(x,y) being 
determined by the specified pressure distribution. Over 
all of the xy plane, Au 0 = 0 applies. 

In all cases, induced velocity u is related to pressure coef- 
ficient Cp by the relation 



Unsteady state .— The steady— state boundary conditions have been 
given in the most general terms possible. The unsteady— state condi- 
tions will be limited to a more restricted type of problem, namely, 
cases wherein the airfoil is assumed to experience at t=0 either an 
abrupt change in angle of attack without pitching or starts to travel 
at the instantaneous velocity V 0 and angle of attack a. In this 
way the transient variation of load distribution and airfoil charac- 
teristics can be calculated for unit angle— of— attack change. Similar 
methods can treat unit rate of pitching, or deflection of aileron, as 
well as the effects produced 
when the airfoil enters a gust 
of given structure . The use 
of solutions of such problems 
in connection with operational 
methods is well known in applied 
mathematics. Applications of 
these operational methods to 
aerodynamics have been given by 
R. T. Jones (references 5 and 6) 
for incompressible fluid theory 
and by Heaslet and Lomax (refer- 
ence 7) for supersonic flow. 

If the rectangular coor- 
dinate system x^z^t* assoc- 
iated with equation (2) is con- 
sidered to be fixed, the airfoil 
moves in the negative x* direc- 
tion and the free— stream velocity 
is zero. A simple distortion of 
the time axis is also introduced 
to simplify the differential Figure 1.— Diagram for use in deter- 

equation. Figures 1(a) and 1(b) mining boundary conditions in two- 

aid in the visualization of the dimensional unsteady motion. 



(a) Supersonic wing. 



(b) Subsonic wing. 
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problems Involved. The airfoil section is assumed to lie initially 
on the x f axis with leading edge at the origin and trailing edge at 
the point x* = c 0 . As time progresses the airfoil sweeps across a 
portion of the x’t’ plane, the leading edge traversing the line 
x* = — Mot’ while the trace of the trailing edge is the line 
x* = c 0 — Mot* . The region hounded hy these lines and the line 
t’ = 0 is that swept hy the airfoil. The characteristic cones of 
the differential equation cut the x’t’ plane along lines inclined 
at ± 45° to either axis . If the airfoil experiences an angle— of— 
attack change a without pitching, the "area" swept over hy the wing 
must yield w = — V Q a. On the other hand, if the airfoil enters a 
gust of constant vertical velocity w 0 , the region over which the 
modification of w is effective is restricted to the region occupied 
simultaneously hy the airfoil and the gust. If, for example, the 
edge of the gust is fixed along the t* axis, this axis will form the 
right-hand "boundary of the region over which the change in boundary 
conditions occurs. A statement of these boundary conditions may be 
put in the following form: 

1. Lifting surface undergoing abrupt change or starting from 

rest with given velocity.— The condition Auq = 0 holds 
over all of the x’t’ plane except for the region swept 
across hy the airfoil. In this latter region, the rela- 
tion w = f(x’,t’) is given, the function f(x*,t’) being 
determined hy the modification in airfoil angle of attack, 
pitching velocity, aileron deflection, or by the gust 
structure. Over all of the x’t’ plane, Aw 0 = 0 applies. 

The expression for pressure coefficient is 


General treatment .— The use of Green’s theorem in the solution 
of second— order partial differential equations leads one to the con- 
sideration of certain particular solutions of the given equations. 


C. 


T? 



-2 dcp 
VqMq dt* 


Solution of Boundary— Value Problems 


NACA TN No. Ij6j 


9 


Because of the physical importance as well as the mathematical appli- 
cability, attention has been centered on the use of a so-called 
fundamental solution or source potential. Thus, in the subsonic 
case, the potential at the point x, z of a unit source located at 
the point 4, £ and applicable to equation (A), table I, is the 
logarithmic function 

cp(x, z) = ^- In r. = ^- In V(x-|) 2 +P 2 (z-^) 2 (3) 

2tt 2it 


while for equation (C) the potential at x,y, z of the ■unit source 
at is 


9(x,y,z) 


-1 = -1 
4itr c 4it */(x-£ ) 2 +P 2 [(y-T|) +(z-C) 2 ] 


(4) 


Here and elsewhere we have (3 2 = | 1 -Mq 2 | where the bars indicate 
that absolute values are to be taken. 

The application of these potential functions to the solution of 
boundary—' value problems in subsonic linearized flow is well known. 
Supersonic theory, however, introduces added complications when 
the fundamental solutions are considered and, although methods have 
been established, the mathematical techniques are of comparatively 
recent origin. The principal difficulty lies in the integration of 
higher— ordered singularities that appear in the three-dimensional 
analysis. Hadamard (reference 4) resolved these difficulties and 
thus avoided the more specialized approach of Volterra (reference 8). 
It would appear, however, that a more direct method of derivation 
stems from Marcel Riesz’s use of fractional integration. (See, in 
this connection, references 9 and 10.) The oddness or evenness of 
the number of dimensions still involves considerable differences 
but the final solutions are easily applied. 

In two-dimensional supersonic flow, the potential at the point 
x, z of a unit source located at the point |, 4 is defined as 
follows : 


cp(x,z) = 0 for (x-|) 2 < P 2 (z— l;) 2 
cp(x, z) = “t for (x-4) 2 ^ p"(z -4) 2 


( 5 ) 
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In three dimensions, the source potential is 


<P(x,y,z) 


-1 _ -a 

2«r D ~ 2«/(x-0 2 HB 2 [(y-^) 2 +(z-U 2 ] 


(6) 


at all points for which the radical is real and is zero elsewhere. 

These functions are directly applicable to equations (B) and (D) 
of table I. Equations (E) and (F) are, of course, special mathe- 
matical cases of equations (A) and (D) for which Mq is 0 and//2, 

respectively. 

By means of the various source potentials, it is now possible 
to present solutions of the differential equations in terms of the 
prescribed boundary conditions. These conditions are assumed to be 
given in the z=0 plane and subscripts u and l shall again 
denote conditions at z=0+ and z=0— , respectively. The general 
solutions appear as follows: 

Equation (A), Mq< 1, 


" (X ’ Z,= SF £ 


Equation (B), 


in < r A>, 

A £=o \oz oz 


— (cpu.-'PZ ) 


| Znr A 


L )s=0 . 


(7) 


M 0 > 1 


cp(x,z) 



— 00 


for z > 0 


:+3z 


l 

~Sx 


d^ for z< 0 


(8) 


Equation ( C ) , Mq < 1 
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Equation (D) , Mq>1 


cp(x,y,z) 


rr 

r m , 

u 

Scpz\ 

JJ 

fi 

ITT 

II 

O 

^ dz 

Sz / 


(<Pu-9l) (-rr £-) ~j 

\o£ r D/t=oJ 


(10) 


In the last equation, the range of integration is confined to that 
portion of the z=0 plane that lies within the Mach forecone of the 
point x,y,z, that is, within the half-portion of the right circular 
cone 


(x-i) 2 - 3 2 [ (y- t)) 2 +(z-C) 2 ] = 0 


lying upstream of the point x,y, z. The semi vertex angle p of this 
cone is the Mach angle and is given hy the relation 

p = arc sin ~ = arc cot B 

M 0 

The symbol [ was introduced by Hadamard and denotes the 

"finite part" of the integral. As in the case of Cauchy 1 s principal 
value, an improper integral is reduced by a prescribed technique to 
a finite and unique value. By definition (see also reference 11), 


r*° A(x)dx _ r x ° A(x)-A (x n ) _ 2A(x 0 ) 

J a (x 0 -x) 3 '' 2 7 a (x 0 -x) 3 ' 2 (x^) 1 ' 2 


In the two-dimensional supersonic case, the solution for the 
velocity potential is expressed as the integral of a distribution of 
source potentials. In all other cases, both the source potential 
and its derivative appear in the integrand, this latter expression 
being identified with the doublet potential. 

Nonlifting case (symmetrical wing, boundary— value problem of the 
first t kind .- Equations (7), (8), (9)7 and (10) are applicable directly 
to the calculation of the potential function corresponding to a 
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symmetrical wing. The relation 9 U = Tz follows from the condition 

u u - vli . Moreover, if — — denotes the local slope of the upper 

dx 

surface of the wing, 

_L = _ i_ ^Z _ ^£u\ 

V 0 8 z V 0 8 z \ dx / 

and the solutions of the various equations are expressed in terms of 
source distributions alone. 

For example, equation (10) becomes 


cp(x,y,z) 



dz u d£dT) 

^ -/ ( x— £ ) 2 — 3 2 [ ( y— q ) 2 +z 2 ] 


(12) 


where the finite part sign is dropped since the integral is proper. 
This equation was given by Puckett in reference 3. 

The pressure coefficient on the surface of the wing is 


c i> * ■ % (sXo 

Lifting case (boundary—' value problem of the first kind) .— From 
the condition wq = w j we have Aw 0 = 0 and the integrands in 
equations (7)> ( 9 ), and (10) are expressed solely in terms of 
doublet distributions, while equation (8) yields the result that 
conditions on either side of the wing have no effect on the other 
side. 


Taking equation (10) as an example, the solution under the pre- 
scribed conditions is 


cp(x,y,z) = - — 


r r m - 7 — — t— 

Jt J |( X -?) 2 -3 2 [(y-r,) 2 +z 2 ]| 3 / 2 


( 13 ) 
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where 





Au 0 ( £,q)d| 


A more direct evaluation of perturbation velocity u can he obtained 
from the alternate expression 


u 


_ 1 _ 

2 * 



Au Q U,n)|3 2 2 dg dr) 
|(x-£) 2 -0 2 [(y-T)) 2 +z 2 ]| / 


Similar expressions exist for equations (7) and (9). 


( 14 ) 


Lifting case ( boundary— val ue problem of the second kind) This 
type of boundary condition cannot be solved directly by means of the 
formulas which have been presented but resolves always into the 
required solution of an integral equation. In three-dimensional sub- 
sonic wing theory, the method of solution depends usually on same 
modification of Prandtl»s lifting-line theory although, more recently, 
lifting— surface theories by Falkner (reference 12) and Cohen (refer- 
ence 13) have been applied successfully. 


In the case of three-dimensional supersonic— wing theory, sources, 
sinks, and doublets have been utilized in two ways in the solution 
of lifting— surface problems. The first of these methods was given by 
Eward (references lk and 15) and is particularly powerful when one of 
the leading edges of the wing is of the supersonic type, that is, 
when the velocity component of the free stream normal to the edge is 
greater than the speed of sound. A second method of solution was 
presented in reference 16 for the important case of wings with sub- 
sonic leading edges, provided the flow field about the wing is of the 
conical type introduced by Busemann (reference 17). The essential 
feature of this method is the use of differential lifting elements 
carrying a constant load and designed for use in conical flow fields. 
The solution consists of determining the distribution of loading 
over these elements so that the resultant induced vertical velocity 
at any point on the lifting surface satisfies the local boundary 
condition. When approached from this standpoint, the problem again 
requires the solution of an integral equation but the equation is of 
the form 


w(x) 



f(xi)dxi 
xi— x 


(15) 
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and is thus veil known from low-speed airfoil— section theory. 
Inversions of this equation have been provided by Allen (refer- 
ence 18) and von Mises and 
Friedrichs (reference 19) • 

Figure 2 shows the 
elemental lifting surface to be 
used. The sides of the element 
extend back from the tip of the 
Mach cone, making angles S and 
&4AS with the positive x axis 
or free— stream direction. The 
vertical velocity induced at the 
point x,y,0 by the element will 
be a function of 5, AS, and ^ 

or, changing the notation, 0, A0 , 



Figure 2.— Lifting surface element 
carrying constant load. 

and to where 


0 = 3 tan S 
04A0 = (3 tan (S + AS) 

--"I 


Denoting the gradient of this vertical velocity with respect to 0 

^W n —A 

by , it follows that 

80 

dwz=0 _ w(0+A0,co)— w(0,co) 

d 0 A0 ^ 0 A0 

where w(0,co) and w(04A0,co) are the velocities induced by rights 
triangular lifting surfaces with constant loading and with vertex 
angles equal to 5 and &4AS, respectively. The velocities 
induced by the constantly loaded surface are determined directly 
from equation (13)- The results of these calculations yield the 
following expressions: 


For co < 0 
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££_ 7 “ -^ 02 - ** ( 16 a) 

de 2rtV 0 J_ 1 <JQi(a>i— 0) 2 

and for cd > 0 


^5=2 ss_ r m tol da.) 

d0 2 jtV 0 Ji ODi(a>i-0) 2 

where 

2 q 

The term — is the load coefficient and is equal to the difference 

q 

between pressure coefficients on the lower and upper surface of the 
wing, 


4P 

q 


PZ — Pu 


= C PZ “ C P U 


Nonl ifting case (boundary— value problem of the second kind) . — In 
the previously discussed lifting case, the induced vertical— velocity 
field of a constantly loaded element was calculated. An analogous 
type of element can also be developed for use in the determination of 
nonlifting wings with prescribed pressure distributions. It is 
apparent that differential expressions similar to equations (l6a) and 
(l6b) must be derived which establish the induced field of the 
x component of perturbation velocity for a conical— flow element with 
constant vertical velocity. Prom equation (12), the following 
expressions result: 

For a) < 0 


dCp _ 2X r°* coi<3n)i 

d0 PrtJ-! (cdi- 0) 2 Vl-ooi 2 


(17a) 
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and for > 0 


dCp 2X ft midjDx 

<30 pit J (oil— 0) 2 7l^i 2 


(1Tb) 


vhere pressure coefficient Cp and surface slope X are 


Cp 


— X = 

Vo’ V 0 


The application of equations (17a) and (17b) to the determination 
of a thickness distribution supporting a given pressure distribution 
consists of determining X as a function of 0 such that the desired 
pressure results. The essential simplification of the method is brought 
about by the use of elements which lead to single integral equations of 
standard form. 


APPLICATIONS 

Yawed Triangular Lifting Surface 

Consider a yawed, triangular flat plate with subsonic leading 
edges such as is indicated in figure 3* Relative to the x axis or 
free— stream direction, the sides of the triangle make the angles 

5 0 and 6i so that the total 
vertex angle is 5 0 + bi=^. 

The quantities 0 O , Q x , and co 
are also introduced where 


Figure 3*— Yawed triangular flat- 
plate lifting surface with sub- 
sonic leading edges. 


0 O = P tan 5 q 


0i = p tan 6 X 


o> = p^ 

X 


A r A 

N. \ 

'"X \ X 

N I t 
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The vertical induced velocity at any point on the wing can be 
found from equations (l6a) and (l6b), provided the distribution of 
the loading factor C is known. Setting C=C(0), the downwash is 
given by the expression: 


v z=0 = - 


3 

2irV 0 


f C(0)d0 P -V 1 "^! 2 

J— G q J ± ^1(^1“ 9) 2 


da>i 



's/l—ad 1 2 

COi ( COjl — 0 ) ^ 



( 18 ) 


Since the lifting surface is flat, the function C (0) must be found 
such that w z =o is independent of to for — 0i< co< 0 O . The integral 
equation can be greatly simplified by integrating the coi variation 
by parts and then taking the partial derivative of both sides of the 
equation with respect to co. In this manner, equation (l8) reduces 
to 


Q Q 

0 = JL P ° C(0)d0 + 1 r ° C(0)d0 

“ J-e 1 oy-e co 2 Vi-co 2 J_ Qi co-0 

vhich becomes 

0 = A p 6 ° c(e)d 0 

8c 0 J CO-0 

-01 


The solution of equation (20) can be written 




A + 


I 

1=1 


Bi 

(0+Xi) 


J (0+0i) (0— 0 q) 


-01 < \ < 6 0 


(19) 


( 20 ) 


( 21 ) 


and, if the integrated loading of the wing is finite 


C(0) 


A0+B 

a/( 0+0i ) ( 0— 0q) 


( 22 ) 
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Substitution of equation (22) into equation (18) yields the two 
relations 


w z=0 = 27 ^ [AHi(©o^®i) + BH 2(0 o^ 9 i) ] 

and 


_ JL 

Z =0 2 V n 


^<7 —I 


[— AH 1 (0 1 ,0 o ) + BHaCQi^o)] 


where 


Hi(0i,0o) = 


-I 


dcgi 


0)1 J (l— COi 2 ) (a)i— 01 ) (o)i+ 0 o) 


and 


H 2 (en,0o)= 7 , ^ 

J 0>1 2 </ (l— (JOjl 2 ) ( CDi — 0 1 ) (o>L+ 0 o) 


(23) 

( 24 ) 

(25) 


(26) 


Equations (25) and (26) may be integrated by the standard methods 
for elliptic integrals and, after substituting into equations (23) 
and (24) and solving for A and B, we have 


A 


V o Wz=0 

pE* 




(27) 


and 


where 


VqVz=0 

(3E» 


20 o 0: 


/ 


2G 


00+^1 


1+0Q01 — a/ ( l— Op 2 -) {1—6 ± 2 ) 


( 28 ) 


(29) 


and E f is the complete elliptic integral of the second kind with 
modulus \! 1— G 2 . 

The load distribution over the wing can now be calculated from 
equations (22), (27), and (28). It follows that 
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&p _ 2C(9) _ _2a ! 2G 

q Vo 2 3E*V e o +0i 

Typical load distributions over a yawed wing are shown in figure 4 for 
3 tan &i = 0.6 and for 3 tan 6q equal to 0, 0.3> 0.6, and 0 . 9 . 


(0q— 0i)0+20q0x 


(30) 



Figure 4.— Angle— of— attack load distribution over yawed triangular 
plan form, 3 tan 63 . = 0 . 6 . 
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Integration of equation (30) over the surface of the triangular 
wing determines for lift coefficient the expression 


CL = 2aa cosA / 2 _ta£_A (31) 

E' J 3 

where A is sideslip angle and 2A is the angle between the leading 
edges. Equation (31) was derived for wings with subsonic leading 
edges and supersonic trailing edge and consequently is valid only for 
cases for which 


[i + A< 90° 

A + A< p 

A - A> 0 (32) 


Downwash Behind Triangular Wing in Supersonic Flow 



The second application will show how doublet distributions may 
be employed in the calculation of downwash in the wake of an unyawed 
triangular wing with subsonic leading edges. The expression for the 

velocity potential will be 
given in general form, but 
in order to avoid detailed 
analysis the value of down- 
wash is determined only 
along the center line in 
the plane of the wing. 


A plan view of the 
wing and wake is shown in 
figure 5. The load distri- 
bution over the wing is 
found from equation ( 30 ) , 
after setting 0o=0i> to 
be 

Ap _ kd 0 2 ajx 

E* 0 P Veo^-rV 2 


( 33 ) 


Figure 5.— Triangular wing and wake 
together with regions used in cal- 
culating downwash. 
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where E* 0 has the modulus a/i— 0 O 2 • Moreover 


cp = / u dx 


so that the jump In potential cpy— cpj on the surface of the wing is 


2aV r 


•Pu— 97 = 

E’oP 


(34) 


and for points in the wake is 


<Pir c P l = ~~ J 0o 2 Co 2 ~P 2 y 2 
E* 0 P 


(35) 


where c Q is the root chord of the wing. Since the wing and its 
wake form a discontinuity surface for the perturbation velocity 
potential and since for all points on this surface 

dcp u _ d? i 

ciz Sz 


it follows from equation (10) that the velocity potential at an 
arbitrary point x,y, z is given by the relation 


<P(x,y,z) 



(36) 


where r is that portion of the wing and wake forward of the Mach 
forecone from the point x,y, z and (<Pu— 1 Pz) is given by equa- 
tions (34) and (35). 


The value of the downwash aft of the wing and along the x axis 
will be calculated from equation ( 36 ), thus 
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is to be determined. In carrying out these calculations, it is 
necessary to consider two segments of the x axis behind the wing: 

Region A extends from the trailing edge to the point where the 
trailing Mach cones from the tips of the wing intersect the x axis 
and thus includes values of x satisfying the inequality 


The final expressions for downwash in the two regions are 
found, after some manipulation, to be 


c 0 ^x< c 0 ( 1+0 o) 


Region B includes values of x for which 


c o (l+0 o ) £ x 


Region A: 


E 2 -(l-k 2 2 )K 2 

jtE’ 0 k 2 




(37) 



(38) 


where 


induced vertical velocity on the wing 


K 


complete elliptic integral of the first kind 


Ki, K 2 complete elliptic integral of the first kind with 
moduli ki, k 2 , respectively 


E 


complete elliptic integral of the second kind 


NACA TN No. 1767 


23 


complete elliptic integral of the second kind with moduli 
ki, k 2 , respectively 

Cq9q 

X— Co 
X— Cp 

c o 0o 

Figure 6 shows the variation of along the x axis for various 

w o 

values of the parameter 0 O =(3 tan 5 0 . The asymptotic values at x=oo 
are also indicated and can he shown to agree with the values of 
downwash at infinity for a wing with the same span load distribution 
in incompressible flow. The discontinuity in downwash at the trailing 
edge is a characteristic property of supersonic— type trailing edges. 


Ei, E 2 

ki 

k 2 



Figure 6.— Variation of downwash on x axis behind a triangular 
wing plotted as a function of distance in chord lengths. 
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Denoting downwash at this point "by wt> Lager strom and Graham 
(reference 20 ) have shown that 


wt = ^*0-00 ^9) 

v 0 E f 0 


A more detailed development of the results shown in figure 6 has 
heen given in reference 21, 


Two-Dimensional Unsteady Lift Problems 

It has already been pointed out that, in the case of unsteady 
motion in a. two-dimensional compressible— flow field, the linearized 
partial— differential equation for the perturbation velocity potential 
can be transformed into the same form that has been considered in 
solving steady-state problems in supersonic wing theory. This 
immediately indicates the possibility that for certain types of 
boundary— value problems in the unsteady case an analogy can be estab- 
lished with three-dimensional lifting— surface problems. 


As an example, let us consider an airfoil that has been flying 
at supersonic speed and then experiences at t t= 0 an abrupt angle— 
of— attack change without pitching. Since the angle— of— attack change 
is assumed to take place at t f =0, it can be assumed that previous 
to this time the induced velocities of the wing are zero and only 
subsequent perturbations are to be calculated. Throughout the swept 
area in the x f t* plane (fig. l) the vertical induced velocity w is 

constant and equal to — VcA. Elsewhere in the z ! =0 plane there is 

ocp 

no discontinuity in the value of pressure, that is, — — 
uous at z 7 =0. 


is contin— 


Suppose now that the area is a wing plan form and that the 
free stream is directed along the t f axis. The characteristic 
cones of the unsteady problem become the Mach cones of the steady- 
state problem, and the Mach number of the free stream is ^/2 since 
the characteristic lines in the figure are inclined k^° to the axes. 


Moreover, the induced vertical velocity is ^ and the perturba— 

, Sz Sep A 

tion velocity in the free— stream direction is ^7* ^ correspond— 

ence can thus be established between the unsteady problem and a 
three-dimensional lifting— sur face problem. 
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As outlined, the boundary— value problem is one of the second 
kind, that is, w is specified on the wing and Au=Aw^=0 off the 
wing. In our particular example, however, the edges of the wing are 
of the supersonic type and no interaction exists between the two 
surfaces of the lifting plate so that pressures on either side can be 
calculated by the methods used for symmetrical nonlifting wings. 

Thus, from equation (12), for z* >0 


<p(t«,x*,z*) 



dt , 1 dx , 1 

*/ (t*— ■ t’l) 2 — (x 1 — X*!) 2 — Z* 2 


(40) 


and for all z* 


< P(t , ,x , ,z t ) = - cp(t* ,x*,-z*) 

The expressions for the indicial load coefficient are as 

follows: 

Region A (between lines x'=t*, t»=0, and x»=CoHMo■t , ) 


Ap 4a 
q Mq 


(4la) 


Region B (between lines x'=-t*, x*=t, and x*=Co-Mot* ) 


Ap _ 4a 
q Vmq 2 -i 


M n x’+t* VM 2-1 /it 

+ 2 f - + arc sin 


— arc cos —21 — — Y 

« x’+Mot* 


itMo 


/it x» \1 

— + arc sin — 

\2 WJ 


(4lb) 


Region C (between lines x’ = — Mot* and x’ = — t 1 ) 


Ap _ 4a 

q J Mq 2 — 1 


(4ic) 
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The growth of 


Ap 


with time, as obtained from equations (4la), 


(4lb), and (4lc), is shown in the portion of figure 7 designed 
"supersonic." At t*=0 the loading jumps to the value ^ and is 
constant along the entire chord. This value persists throughout the 
previously denoted Region A and thus, with advancing time, moves 
rearward along the chord, leaving the trailing edge at t* = - c -2 y . 
Over the forward portion of the chord the familiar Ackeret type of 
steady— state loading becomes effective, spreading back from the 

leading edge and occupying the entire chord length after t* = ,, C °\ . 

c 0 Mo— 1 

Previous to t* = - a transition region between the two types of 

JMq •“ X 

constant loading exists, and subsequent to this time this transition 

region moves aft and leaves the trailing edge at t* = — . 

M 0 — 1 



Figure J Pressure distributions on wing undergoing sudden angle— of— 
attack change at t*=0. 

For purposes of comparison, the growth with time of the angle— of— 
attack indicial load coefficient for subsonic flight is also shown in 
the part of figure 7 entitled "subsonic." Since in this case the 



MCA TN No. 1767 


27 


lifting— surface analogue involves subsonic leading and trailing 
edges, the analysis requires the solution of a boundary— value problem 
of the second kind. The method of Eward (reference 15) was used to 
obtain the results shown. It is to be noted that the expression 

Ap _ 4a 
q Mo 

holds at t*=0 for all values of Mach number. 

2 Mpt 7 


Co 



Figure 8.— Indie ial lift— curve slope for Mach numbers between 

0 and 1.4 shown to time required to travel 12 half— chord lengths. 

Figure 8 shows the variation of the indicial lift function 
CL a (t’) defined by the relation 


CLa(t') 



(42) 
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as a function of Mach number and half— chords s traveled by the 
airfoil. The curve at 1^=0 was first studied by Wagner (refer- 
ence 22) and R. T. Jones (reference 6) . Since the starting value is 

77—, initially be infinite. Immediately afterward, 

Mo 

however, it assumes the value tc and then rises to the asymptotic 
value 2 tt. At a Mach number of 0.4 the starting value of CL rY (t l ) 
is 10 followed by a decrease for the time required to travel approxi- 
mately one— half chord length and finally a steady rise takes place 

to the asymptotic value — = ■■ ■ - At M o =0.8 the behavior is similar. 

a/ 1-M 0 2 

The dashed portions of the curves were determined from the known 
variations of the functions and were not calculated explicitly. The 
asymptotic values of consistent with the Prandtl-Glauert correc- 

tion become so high, however, with increasing Mach number that the 
assumptions of small perturbation theory are undoubtedly invalidated 
near Mo=l for sufficiently large values of s. The initial portions 
of the subsonic curves shown in the figure are, however, valid results 
of the theory. The nature of the indicial lift function is somewhat 
different at supersonic Mach numbers in that the beginning portions of 
the curves are flat. The curves rise afterwards, however, in a finite 
time to their steady— state value. From equations (4la), (4lb), and 
(4lc) the expressions for CL a (t’) are easily calculated for Mo> 1 
and are as follows: 

First time interval 0 < t* < 

1+Mq 

Cl u (t*)=^- (*3a) 

a M q 

Second time interval C ° < t* < C ° 

1+M 0 M 0 -l 


CLjf) = £ 


_ Mo 1,2 + 


arc sin 


Cq-Mo*» 

t» 


J Mq 2 — 1 


arc cos 


t * +M q c Q — t 1 Mq 2 
Co 


+ rr^— ^t* 2 — (co— t’M 0 ) 2 
MqCq 


(^3t) 
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Third time interval 


Co 

Mo-1 


<t* 


CLjt’) 


4 

a/Mq 2 — 1 


(43c) 


Some of the above results, along with further developments 
involving the entrance at supersonic speed of an unrestrained air- 
foil into a gust, have been given in reference 7* 


Swept— Back Wings at M 0 =l 

We turn now to the special form of the basic differential equa- 
tion for the case Mq=1. As given in table I, equation (E), the 
velocity potential satisfies Laplace’s equation 


4 ) yy + cp zz 0 


in two dimensions. The boundary conditions need, therefore, to be 
given along strips normal to the free— stream direction. Equation ( 7 ) 
expresses the solution of the equation in terms of two-dimensional 
sources, sinks, and doublets where now 


Zn r A = Zn V (y- T1 ) 2 +( z -^) 2 


The proposed problem is the determination of the angle— of— 
attack load distribution over a swept— back lifting plate, the 
leading edges will be assumed straight lines while the trailing edge 
will, for the time being, be left arbitrary. The nature of the wing 
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is thus indicated somewhat arbitrarily in figure 9( a ) • 



(b) Plan form satisfying Kutta— Joukowski 
condition. 

Figure 9-— Swept— back wings for analysis at Mq=1. 

Denoting the semivertex angle by S Q so that the equations of the 
leading edges are 


y = ± x tan 5 0 = ± mx 


m 


it follows from equation (7) that since 
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the velocity potential is given hy the relation 


where 


<P(x,y,z) 


z 

2n 



Acpp (x,n) 
(y-Ti) 2 +z 2 


dti 


Acpo(x,y) = 


Av o (x,r),0) dq 


' mx 


( 45 ) 


It is then possible, after integrating equation (45) by parts and 
imposing the condition that A%(x,y) =0 at T) = ± mx, to calcu- 
late the derivative of cp with respect to z and thus obtain for 
Wq, induced vertical velocity on the wing, the expression 




AVq(x,T|) 

y-T] 


dT] 


( 46 ) 


This integral equation is to be solved for Av Q , the velocity w Q 
being assumed constant on the wing and equal to — V Q a where a is 
angle of attack. The load distribution can then be calculated from 
this solution by means of the relation 


Ap _ 2 2 d 

q V 0 8x V 0 dx 


I Av Q (x,q) 

^mx 


dr) 


( 47 ) 


The remainder of the analysis can best be divided into two parts: 
The first case treating values of x between 0 and c Q , the second 
dealing with the remaining values of x on the wing. 

Case I: 0 ^ x < Co.— Since the leading edges of the wing are 
of subsonic type, singularities in pressure occur at these edges so 
that the required solution of equation ( 46 ) is of the form 


A+By 

J m 2 x 2 -y 2 


( 48 ) 


Av 0 (x,y) = 
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Substitution into equation (46) and use of the fact that A%(x,y) is 
an even function of y leads to the result A=0, B=— 2wq. Hence 


Acp 0 (x,y) = - 2w 0 vk^x 2 — y 2 


m 


and 


Ap _ — 4vom 2 x _ 4am 2 x 


q V Q Vm 2 x 2 -y 2 V m 2 x 2 -y 2 


(50) 


Case II: c Q ^ x.- Let the equation of the trailing edge be 


= a(x) or x = a*(y) 


(5D 


Using in equation (46) the fact that 


Av 0 (x,y) = - Av 0 (x,-y) 


the expression for v Q becomes 


-1 r^( x ) 2r}Av 0 dT] l / f' mx 2r>Av 0 dii 

- s j o - a? 7 a(x) 


If, on the surface of the wing, Acp 0 (x,r]) is known, then, in 
the wake, the discontinuity in the velocity potential is AcPo£ a *( T l) j 1 )) 1 
since no contribution to the jump in potential is made past the 
trailing edge. It follows that if on the wing 


Av 0 (x,T)) = f(x,q 2 ) 


(53) 


then, in the wake 
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Av 0 (x,t}) 


A ^Pq(x,ti) 3a/* 
3x 3y 


+ A 


^Po(x,Tl) 

hr] 


x=a *(i]) 


(53) 


Vq Ap [a*(q) ,r|3 3a* 
2 q 3y 


+ f [a*(q) ,r\ 2 ] 


(54) 


Substituting equations (53) and (54) into equation (52) and intro- 
ducing the Kutta condition that loading at the trailing edge is 2 er o ^ 
we get the modified integral equation 



where the variables 
The function 


f[a*( //5i) ,cTi ] 1 o m2x2 f(x,a!) 


0-0 1 


doi- 


2it 


'a 2 


d- hi 


dui 


a, oi now replace y 2 , t) 2 , respectively. 

t(x,a) = a» 0 AgLM 

V m^x^— cr 


(55) 


(56) 


satisfies equation (55) and it remains now to determine a(x) so 
that pressure is zero at the trailing edge. But from equation (53) 


Acp(x,y) = 2w q 


/vi 2 ~ a2 ( x ) 

Jhy v m 2 x 2 — yi 2 


dyi 


and thus 


Ap(x,y) 


= 4a 


da 

— k* — F(k,t) + mE(k,t) + 
dx 


y A^-k’ 2 ] 
x v m 2 x 2 — ; 


— k ,2 m 2 x 2 

•y 2 


(57) 


where k’ = ^ and E(k,\|r) ,F(k,i|r) are incomplete elli ptic integ rals 
with modulus k = *A— k* 2 and argument t=arc sin /l 


m 2 x 2 


At the trailing edge y^a(x) 


and 
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where the elliptic integrals are now complete. If load coefficient 
is set equal to zero, the differential equation 


da _ m 2 x E 
dx a K 


(58) 


where k = / 1 - — , follows . 

v m 2 x 2 

The integration of equation ( 58 ) leads to the shape of the 
trailing edge for which the Kutta condition is satisfied. Figure 9(h) 
shows the plan form of the wing. It can be shown that the slope of 
the extended trailing edge approaches the slope -of the leading edge. 



Figure 10. Load distribution over swept— back plan form at Mq=1. 

In figure 10 the load distribution is shown at three spanwise 
stations for the case when the wing is cut off along a line normal 



NACA TN No. Ij6j 


35 


to the free— stream direction. Over the center section of the wing 
the Ackeret type of distribution exists. The remaining sections have 
discontinuities in slope of the loading at the point where the chord 
line is cut by the Mach cone arising at the trailing edge of the root 
chord. This behavior of the loading has been noted elsewhere for 
swept— back wings at higher Mach numbers. (See, e.g., reference 23.) 

Lift coefficient Cl of the wing is given by the expression 


where . S and b are, respectively, area and semi span of the w in g and 
the first integral extends from the leading edge to the trailing edge. 
This equation may be rewritten as 


where Acp(T.E.,y) is the jump in potential at the trailing edge and 
thus equal to the circulation function r (y) . The following results 
are obtained: 




o 


(59) 


r(y) = 2mV 0 c 0 a, O^y^A 


r(y) = 2V 0 ab E(k,\|r)-(l-k 2 )F(k,\|f) , A^y<~b 


(6o) 


where k = 


/ 



and A is the lateral distance to the inboard 
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tip of the wing. 


(See fig. 9(t>) - ) 



In figure 11 the value of 
r 

is plotted as 

V 0 c 0 a. tan 5 0 

y 

a function of - 

Cq tan o 0 

for a wing with semi span 
b=1.325 c 0 tan 6 0 and aspect 
ratio AR=4.57 tan 5 0 . Results 
of the integration of equation 
(59) a **e shown in figure 12 

CL a 

where - — is plotted as a 

tan S 0 


function of 


AR 

tan S 0 


Figure 11.— Spanwise distribution 
of circulation for swept— back 
wing at Mo=l, b=1.325 c 0 tan S Q 
and AR=1 j-. 57 tan 5 0 . 



AR/TAN 8 0 


Figure 12.— Lift— curve slope as a 
function of aspect ratio for 
swept— back wing at Mq=1. 


The methods presented here 
can be applied to the case of 
the swept— back wing with tips 
cut off parallel to the free 
stream. In this case a Mach 
cone originates not only at the 
trailing edge of the root chord 
but also at the intersection of 
the leading and the lateral 
edges. On the portion of the 
wing downstream of this Mach 
cone, the load distribution is 
modified so that an abrupt dis- 
continuity exists at the Mach 
cone and negative loading is 
effective over this part of 
the wing. A similar effect on 
this type of swept— back wing has 
also been noted at higher Mach 
numbers in reference 23. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


NACA TN No. Ij6j 


37 


a.Q 

A 

AR 

b 


Co 


C, C(0) 

CL 

C La 

CL a (t') 


E 

E», E’q 


Ei, E 2 


E(k,t) 

E(k,t) 

G 


ki 


k 2 


APPENDIX 


List of Important Symbols 
speed- of sound in free stream 

lateral distance to inboard tip of swept wing (See fig. 9) 
aspect ratio 
semi span 

chord length (two dimensions) 
root chord (three dimensions) 

load distribution factor introduced in equation ( l6) 
lift coefficient 


lift— curve slope 


indicial lift function 
pressure coefficient 

complete elliptic integral of the second kind with modulus k 

com plete el liptic integrals of the second kind with moduli 
Vi— G 2 , J 1— 0 O 2 , respectively 

complete elliptic integrals of the second kind with moduli 
ki, k 2 , respectively 



o* , 

/ v 1— k 2 sin 2 cp dcp 

j o 

L 


dcp 

a/i— k 2 sin 2. 


parameter defined in equation (29) 

C qQq 

x- C q 


X— Cn 
c o 0 o 
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K complete elliptic integral of the first kind with modulus k 


Ki, K 2 complete elliptic integrals of the first kind with moduli 
ki> k 2 , respectively 


m 


Mo 

P 

Po 

q 

r A 

r C 


tan 5 0 

free— stream Mach number 
local static pressure 
free— stream static pressure 
free— stream dynamic pressure 
v /(x-i) 2 +p 2 ( z-(;) 2 
v /(x-|) 2 +3 2 [(y--n) 2 +(z-C) 2 ] 



J (x— £ ) 2 — p 2 [ (y— t]) 2 +(z— t,) 2 ] 


s distance traveled in half— chords 

S area of wing 

t time 

u,v,w perturbation velocity components parallel to x,y, z axes, 
respectively 

Au 0 , Awo jump in value of u, w at the z=0 plane 

V Q velocity of free stream 

x,y, z Cartesian coordinates 

x* , z* ,t T coordinates introduced in equation (2) 
a angle of attack in radians 

J 1 1- Mq 2 | 


3 

Ap 

T 


load coefficient 


C P °Pu 
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A 

6 

^o* 




A 

u 


S*n*C 


Po 


T 

<P 


semivertex angle of yawed triangular wing 

angle between lifting element and x axis 

angles between leading edges of yawed triangle and x axis 

0 tan 5, 0 tan 5 0 , 0 tan 6 X 

sideslip angle of yawed triangle 

Mach, angle ( arc sin ) 

\ “o/ 

Cartesian coordinates 
free— stream density 

region of integration in equation ( 10 ) 
perturbation velocity potential 


CD 0^ 

K x 

sign denoting "finite part" of integral 

Subscripts 

u subscript denoting value of variable on upper surface 

of wing 


subscript denoting value of variable on lower surface 
of wing 
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